SCIENTIFIC 

REPORTS 




OPEN 



SUBJECT AREAS: 
SPECTROSCOPY 
OPTICAL METROLOGY 
FREQUENCY COMBS 



Received 
26 November 2013 

Accepted 
2 January 2014 

Published 
22 January 2014 



Correspondence and 
requests for materials 
should be addressed to 
T.Y. (yasui.takeshi® 
tokushima-u.ac.jp) 



Spectrally interleaved, 
comb-mode-resolved spectroscopy using 
swept dual terahertz combs 

Yi-Da Hsieh', Yuki lyonaga 1 , Yoshiyuki Sakaguchi', Shuko Yokoyama 1 ' 2 , Hajime Inaba 3,4 , 
Kaoru Minoshima 4,5 , Francis Hindle 6 , Tsutomu Araki 1 & Takeshi Yasui' A7 

'Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan, 
2 Micro Optics Co., Ltd, 2-20, Oe-Nakayama, Nishikyo, Kyoto 610-1 104, Japan, 3 National Metrology Institute of Japan, 
National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8563, Japan, 
4 ERATO Intelligent Optical Synthesizer Project, JST, 1-5-1 , Chofugaoka, Chofu, Tokyo 1 82-8585, Japan, 5 Graduate School of 
Informatics and Engineering, The University of Electro-Communications, 1-5-1 , Chofugaoka, Chofu, Tokyo 1 82-8585, Japan, 
6 Laboratoire de Physico-Chimie de I'Atmosphere, Universite du Littoral Cote d'Opale, 1 89A Av. Maurice Schumann, Dunkerque 
59140, France, institute of Technology and Science, The University of Tokushima, 2-1 Minami-Josanjima, Tokushima, 
Tokushima 770-8506, Japan. 

Optical frequency combs are innovative tools for broadband spectroscopy because a series of comb modes 
can serve as frequency markers that are traceable to a microwave frequency standard. However, a mode 
distribution that is too discrete limits the spectral sampling interval to the mode frequency spacing even 
though individual mode linewidth is sufficiently narrow. Here, using a combination of a spectral 
interleaving and dual-comb spectroscopy in the terahertz (THz) region, we achieved a spectral sampling 
interval equal to the mode linewidth rather than the mode spacing. The spectrally interleaved THz comb was 
realized by sweeping the laser repetition frequency and interleaving additional frequency marks. In 
low-pressure gas spectroscopy, we achieved an improved spectral sampling density of 2.5 MHz and 
enhanced spectral accuracy of 8.39 X 10" 7 in the THz region. The proposed method is a powerful tool for 
simultaneously achieving high resolution, high accuracy, and broad spectral coverage in THz spectroscopy. 

Use of optical frequency combs has emerged as a new technique for broadband spectroscopy because they 
have the potential to achieve high spectral resolution in a broad spectrum in addition to high accuracy 
frequency markers 13 . Optical comb has been usually used in combination with a Michelson-interfero- 
meter-based Fourier-transform spectrometer (FTS) 4 . However, the spectral resolution in usual FTSs is insuf- 
ficient to spectrally resolve each mode because the modes are distributed too densely for FTS. 

Recently, dual-comb spectroscopy has unlocked further potential for optical-comb-based broadband spectro- 
scopy 512 . In this method, one optical comb that contains the measurement information (frequency interval = 
frepi) is rea d with another optical comb with a slight frequency offset (frequency interval = f rep2 = f repl + /offset)- 
This results in the generation of radio-frequency (RF) comb (frequency interval = /offset)' Since the RF comb is 
essentially a replica of the measurement comb, downscaled in frequency by the ratio / repl // 0 ff set , one can utilize the 
comb-mode-resolved spectrum via direct observation of the RF comb. Since the first experimental demonstration 
of the dual-comb spectroscopy in mid- infrared region 6 and the following observation of the comb-mode-resolved 
spectrum in terahertz (THz) region 7 , the implementation of dual-comb spectroscopy has made notable progress 
in various wavelength regions because of high sensitivity, rapid data acquisition, high spectral resolution, and 
excellent accuracy 8-12 . 

The optical comb mode enables us to achieve a frequency resolution equal to the comb mode linewidth 
(typically, less than 1 MHz) when the absolute frequency is determined based on a single comb mode. 
However, when the comb-mode-resolved spectrum obtained by dual-comb spectroscopy is used for broadband 
spectroscopy, the spectral sampling interval is limited to the comb mode spacing (typically, from several tens to a 
few hundred MHz) due to the excessively discrete mode distribution. If the comb mode gap is filled by inter- 
leaving additional frequency marks, the spectral sampling density and accuracy would be vastly enhanced in the 
broadband spectroscopy. The combination of such interleaving with dual-comb spectroscopy has been demon- 
strated in optical region 1314 and sub-Doppler spectroscopy has now been achieved by making full use of resolution 
improvement through interleaving 15,16 . However, there have been no attempts to interleave THz comb. 
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In this article, we propose a spectrally interleaved dual-THz- 
comb spectroscopy based on a swept dual-THz-comb technique. 
Furthermore, the proposed method was effectively applied to high- 
resolution, high-accuracy, broadband spectroscopy of gas-phase 
molecules at low pressure. 

Results 

Principle of operation. We first describe spectral behavior of dual- 
THz-comb spectroscopy in optical, THz, and RF regions as shown in 
Fig. 1(a) 7 . When the first optical comb (frequency interval = f repl ) is 
incident onto an electrically-biased photoconductive antenna for 
THz emission (PC A emitter), it is down-converted to the THz 
region without any change to the frequency spacing because the 
PCA emitter squares the electric field of the optical comb via the 
quadratic detection, converts the optical intensity to electric carriers, 
and emits a frequency comb of electromagnetic THz wave (EM-THz 
comb) via the dipole radiation. The EM-THz comb is a harmonic 
frequency comb of f repl without a frequency offset (freq. = f rep i, 
2frepi> 3f repl , nf repl ). Next, we consider what happens when 

the second optical comb (frequency interval = f rep2 = f rep i + /offset) 
slightly detuned from the first one is incident onto a non-biased, 
photoconductive antenna for THz detection (PCA detector). In 
this case, THz comb of photocarrier with frequency spacing oif rep2 



[PC-THz comb; freq. = f rep2 , 2f rep2 , 3f rep2 

/offset)* 2(f repl + foffset)> 3(f repl + f 0 ff s et)> " 



n frep2 = (f, 



repl 



'•*•. "(/repi + /offset)] IS 

induced in PCA detector. Detection of EM-THz comb with the 
PCA detector including the PC-THz comb allows us to detect the 
RF comb (frequency interval = f 0 ff S et) as a result of the multi- 
frequency-heterodyning photoconductive mixing occurring 
between the EM-THz and PC-THz combs. 

Frequency-domain signal acquisition using an RF spectrum ana- 
lyzer enables a direct observation of comb-mode-resolved spectrum 
in RF comb 7 . However, signal acquisition efficiency of RF comb is 
limited by super-heterodyne technique used in RF spectrum analyzer 
because it acquires only the single frequency-component signal coin- 
ciding with the frequency of its local oscillator at any given moment 
and then obtains the overall spectrum by sweeping the frequency of 
the local oscillator. On the other hand, a combination of time- 
domain signal acquisition with Fourier transformation (FT) 
improves signal acquisition efficiency because it can capture all fre- 
quency components of the RF comb at once in the time domain. To 
obtain the THz-comb-mode-resolved spectrum at better signal-to- 
noise ratio and dynamic range, we used an asynchronous-optical- 
sampling THz time-domain spectroscopy (ASOPS-THz-TDS) 



system consisting of two mode-locked lasers with slightly mis- 
matched repetition frequencies 11,17-20 . Although the ASOPS-THz- 
TDS has been usually used for acquiring a temporal waveform of a 
single THz pulse within one pulse period (= l/f rep i), the non-mech- 
anical nature of the time-delay scanning enables us to extend the time 
window t w up to multiple pulse periods (=N p /f repl , where N p is 
number of THz pulse) in order to acquire the temporal waveform 
of a THz pulse train, as shown in Fig. 1 (b) . Taking FT of the temporal 
waveform gives the detailed spectrum of the THz comb modes 
because FT of the periodical THz pulses imprints frequency modu- 
lation on the broadband THz spectrum. We call it a standard THz 
comb. In this case, frequency spacing between comb modes is equal 
t° frepi whereas linewidth of the comb mode is determined by l/t w , 
namely, f repl /N p . 

We should consider how we could interleave the additional marks 
into the gap produced between the comb modes of a standard THz 
comb. The THz comb is a harmonic frequency comb of the laser 
repetition frequency f repl and does not contain a carrier-envelope- 
offset frequency f ceo . Therefore, the absolute frequency of each comb 
mode ( = mf repl , where m is the order of THz comb modes) can be 
tuned by changing/ repJ . For example, since a single sweep off repl by 
dfrepi causes a frequency shift of m-order mode (freq. = mf repl ) by 
m3f repl , M-times repetition of this sweep will lead to a total shift of 
MmSf repl in mf repl . If incremental sweeping of the comb mode is 
repeated at an interval equal to the mode linewidth, namely, mSf repl 
= frepi/N p or Sf repl = f re pi/(tnN p ), and all of the resulting comb 
spectra are overlaid in the spectral domain, as shown in Fig. 1(c), 
the frequency gaps of a standard comb can be completely removed. 
To fill the frequency gaps, such incremental sweeping has to be 
repeated times of M l=f r e P i/(mSf repl ) = N p ]. In this way, a spectrally 
interleaved THz comb will be achieved. This is equivalent to con- 
tinuous sweeping of a single-mode, narrow-linewidth CW-THz 
wave. The resulting spectral resolution will be equal to the linewidth 
of the comb mode. It is important to note that M is sufficiently large 
to fully interleave the comb gap because the tunable range oif repl is 
within 1% of/ repJ at maximum. 

Figure 2 illustrates a schematic diagram of the experimental setup, 
which contains dual mode-locked Er-doped fibre lasers (ASOPS 
TWIN 250 with P250, Menlo Systems; centre wavelength a c = 
1550 nm, pulse duration At = 50 fs, mean power P mean = 
500 mW,/ repJ = 250,000,000 Hz,/ rep2 = 250,000,050 Hz, and/ 0#se( 
= 50 Hz) and a THz optical setup for low-pressure gas spectro- 
scopy 12 . The temporal waveform of the pulse train of the THz electric 
field was acquired at a sampling interval of 100 fs. Finally, the 
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Spectral behavior of dual-THz-comb spectroscopy in optical, THz, and RF regions, (b) Temporal waveform of THz pulse train and 
I THz comb spectrum, (c) Spectrally interleaved THz comb achieved by incremental sweeping of THz comb mode and spectral 
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Figure 2 | Experimental setup. Rb-FS, rubidium frequency standard; SFG-XC, sum-frequency-generation cross-correlator; SHG, second-harmonic- 
generation crystals; L, lenses; EDFA, erbium-doped fibre amplifier; OSC, erbium-doped fibre oscillator; PCA1, dipole-shaped low-temperature-GaAs 
photoconductive antenna for THz emitter; PCA2, dipole-shaped low-temperature-GaAs photoconductive antenna for THz detector; Si-L, silicon lenses; 
AMP, current preamplifier. 



amplitude spectrum of the THz comb was obtained by FT of the 
temporal waveform of the THz pulse train. 

Spectroscopy of water vapour. Since water vapour plays an impor- 
tant role in the Earth system, particularly in the radiation budget, it is 
important to measure it correctly from the viewpoint of global 
warming 21 . Here, the rotational transition 1 10 <— loi at 0.557 THz 
in water vapour was measured to assess the system capacity to resolve 
fine spectral signatures 22 . A sample gas, a mixture of water vapour 
(partial pressure =10 Pa) and nitrogen (partial pressure = 320 Pa), 
was introduced into the low-pressure gas cell. The sample has an 
expected pressure-broadening linewidth of 23 MHz full-width at 
half-maximum (FWHM) from the self-broadening of water 
vapour, the collision-broadening induced by nitrogen 21 , and an 
expected Doppler-broadening linewidth of 0.026 MHz 23 . 

Figure 3(a) shows a power spectrum of the whole standard THz 
comb, obtained by taking the FT of the temporal waveform of 10 
consecutive THz pulses {N p = 10, t w = Np/f re p; = 40 ns, number of 
signals averaged = 5000, acquisition time = 1000 s). The maximum 
dynamic range of 40 dB in power was achieved at 0.5 THz, and the 
observed spectral bandwidth was extended up to 2 THz. Figure 3(b) 
shows the amplitude spectrum of the standard THz comb expanded 
around 0.557 THz. The comb modes had a frequency gap of 
250 MHz and a linewidth of 25 MHz. This mode gap was exactly 
equal to f repl as shown in Fig. 1(b). On the other hand, the mode 
linewidth here was determined by the reciprocal of the temporal 
window t w (=40 ns), although the intrinsic linewidth is as little as 
a few Hz 24 . The amplitude spectrum without incremental sweeping 
of the THz comb modes did not indicate the spectral shape of the 
absorption line due to the excessively coarse distribution of the comb 
modes compared with the absorption linewidth of the low-pressure 
water vapour. 

Next, we demonstrated incremental sweeping of the comb modes 
across the absorption line at 0.557 THz by changing the repetition 
frequencies of both lasers. Incremental increases of f repl and f rep2 
{=fre P i + /offset) by 11,220.8 Hz (=Sf repl = 5f rep2 ) were repeated 
ten times (M = N p = 10) while keeping f a g set at 50 Hz. The resulting 
overlaid spectra of the comb modes are shown in Fig. 3(c). In this 
demonstration, a single shift oif repl by 0.004488% resulted in sweep- 
ing of the comb modes by 10% of their interval ( = 25 MHz) due to 
the large number of THz comb modes (m = 2,228 at 0.557 THz). 



The frequency gaps between the comb modes in Fig. 3(b) were filled 
by interleaving additional frequency marks. In other words, a spec- 
trally interleaved THz comb was successfully achieved. Although the 
comb mode amplitude is undulated due to a slight mismatching 
between the comb mode peaks and the spectral sampling points, a 
sharp spectral dip clearly appeared at the position of the water 
absorption line of 1 10 <— 1 01 . The absorption was extracted from 
the power spectrum by taking the peak amplitude of each mode 
normalized with a reference spectrum obtained under identical con- 
ditions, as shown in Fig. 3(d). To determine the spectral linewidth, we 
fitted a Lorentzian function to the measured spectral profile because 
the spectral shape in THz region arises from the pressure broadening 
rather than Doppler broadening ( = 0.026 MHz) 23 at the present 
experimental condition. As a result, the spectral linewidth was deter- 
mined to be 24 MHz, indicated by a red solid line in Fig. 3(d), and 
this value is consistent with the expected pressure broadening line- 
width ( = 23 MHz). 

The same procedure was used to examine the 1 10 <— loi water 
absorption line for various pressures. The black circles in Fig. 3(e) 
show FWHM of the observed absorption line with respect to the 
partial pressure of the water vapour, which was varied between 
5 Pa and 160 Pa. For comparison, the pressure broadening char- 
acteristic of this gas sample is expected from the self-broadening of 
water vapour and the collision broadening induced by nitrogen, as a 
red solid line in Fig. 3(e) 21 . The experimental data were in good 
agreement with the expected line down to 25 MHz, and then 
deviated from it. This result clearly indicated that the increased spec- 
tral sampling density in the interleaved THz comb enhanced the 
spectral resolution of this gas spectroscopy to the level of the line- 
width of the comb mode ( = 25 MHz) rather than the gap size 
( = 250 MHz) for unswept measurements. 

Spectroscopy of gas-phase acetonitrile. Since acetonitrile (CH 3 CN) 
is a very abundant species in the interstellar medium and is an ideal 
observational probe of the kinetic temperature and density of 
interstellar clouds 25 , we performed gas-phase spectroscopy of this 
molecule. Since CH 3 CN is a symmetric top molecule with a 
rotational constant B of 9.2 GHz, it displays two features: a series 
of manifolds of rotational transitions regularly spaced by 2B 
( = 18.4 GHz), and hyperfine structure of rotational transitions into 
each manifold determined by the centrifugal distortion constant 
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Figure 3 | (a) Power spectrum of the whole standard THz comb. Amplitude spectra of (b) standard THz comb and (c) interleaved THz comb around 
0.557 THz after passing through low-pressure water vapour contained in the gas cell, (d) Absorption spectrum and (e) pressure broadening 
characteristics of the rotational transition 1 10 «— 1 01 in water vapour. 



Dj K 25 . Conventionally, it has been necessary to measure these two 
features separately by using broadband THz-TDS 26 and a high- 
resolution CW-THz spectrometer 27 due to the limited frequency 
resolution, accuracy, and/or coverage in traditional THz 
spectroscopy. Recently, one-pulse-period ASOPS-THz-TDS was 
used to observe both the manifold and the hyperfine-structure 
features simultaneously 20 . However, since its spectral resolution 
remained at the laser repetition frequency, it was still difficult to 
fully resolve all of the hyperfine structure. The reduced comb- 
mode linewidth and increased spectral sampling density offered by 
the spectrally interleaved dual-THz-comb configuration should 
allow the hyperfine structure to be resolved at high resolution and 
accuracy. 

First, we performed spectroscopy of gas-phase CH 3 CN at 40 Pa 
enclosed in a gas cell by using the standard THz comb, obtained by 
taking the FT of the temporal waveform of 10 consecutive THz pulses 
(N p = 10, t w = N p /f repl = 40 ns). Figure 4(a) shows the absorption 
spectrum of this gas sample obtained by using the standard THz 
comb, in which a series of manifolds from J = 16 to J = 55 was 
clearly confirmed at intervals of 18.4 GHz within a frequency range 



from 0.3 to 1.0 THz in the same manner as the THz-TDS 26 . Next, we 
used an interleaved THz comb with a mode linewidth of 25 MHz for 
high-precision spectroscopy and allowed the hyperfine structure into 
a single manifold around 0.64 THz (J = 35 to J = 34) to be observed, 
as shown in Fig. 4(b), in the same manner as the high-resolution CW- 
THz spectrometer 27 . To assign these absorption lines, multi-peak 
fitting analysis based on a Lorentzian function was used, indicated 
by the red solid line in Fig. 4(b). In comparison with literature values 
reported in the JPL database (see blue dashed line) 28 , we successfully 
assigned lines K = 2 to 10 within a frequency discrepancy of 4.7 ± 
2.6 MHz (mean ± standard deviation for 9 absorption lines). To 
assign lines K = 0 and 1 with the frequency separation of 12 MHz, 
the linewidth and incremental step of the comb modes should be 
further reduced. To this end, we expanded the time window t w to 
400 ns (=Np/f rl ,p 1 ), in which the temporal waveform of 100 consec- 
utive THz pulses was measured (Np = 100). It should be emphasized 
that the time window of 400 ns corresponds to a mechanical time- 
delay scanning by 60 m when a traditional THz-TDS is used. Only 
ASOPS-THz-TDS can provide such an incomparably wide time win- 
dow. This reduces the linewidth of the comb mode to 2.5 MHz. 
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Figure 4 | Absorption spectrum of low-pressure acetonitrile gas (a) within a frequency range from 0.3 to 1 THz and (b) around 0.6428 THz. (c) 
Amplitude spectrum and (d) absorption spectrum obtained by 15 incremental sweeps of THz comb mode across two adjacent absorption peaks. 



Incremental sweeping of this comb mode was repeated 15 times at 
intervals of 2.5 MHz across the spectral region of the two absorption 
lines of K = 0 and 1 (gas pressure = 20 Pa, pressure broadening 
linewidth = 5.9 MHz, Doppler broadening linewidth = 
0.020 MHz) 24,29 . Figure 4(c) indicated that two spectral dips clearly- 
appeared although the pressure broadening makes them overlap 
partially. The spectral linewidth was determined to be 7.5 MHz for 
K = 0 and 6.1 MHz for K = 1 by performing multi-peak fitting 
analysis, as shown in Fig. 4(d), which are consistent with the expected 
pressure broadening linewidth ( = 5.9 MHz). These results clearly 
indicated that the reduced comb-mode linewidth and the increased 
spectral sampling density further enhanced the spectral resolution. 
We also determined their centre frequencies to be 0.643257578 THz 
for K = 0 and 0.643269502 THz for K = 1 from the multi-peak 
fitting analysis. The discrepancy of the centre frequencies from the 
literature values in the JPL database (see blue dashed line) 28 was 
0.578 MHz for K = 0 and 0.502 MHz for K = 1, corresponding to 
a mean spectral accuracy of 8.39 X 10~ 7 . Since these discrepancies 
were within the spectral sampling interval of 2.5 MHz, further reduc- 
tion of the mode linewidth and incremental step will improve the 
spectral accuracy as well as the spectral resolution. These results 
clearly indicated the high potential of the spectrally interleaved 
dual-THz-comb spectroscopy to simultaneously probe multiple 
absorption lines of low-pressure molecular gas. 

The accuracy of the absorption strength is also important to make 
a reliable quantitative analysis in gas-phase spectroscopy under low 
pressure. To evaluate the quality of the obtained spectrum, we made a 
simple comparison of our spectrum with the integrated absorption 
intensities in JPL database 28 , indicated as blue solid lines marked with 
blue cross in Fig. 4(b). The integrated absorption intensity is defined 
as an integration of the absorption spectrum and depends on both 
the absorption coefficient and linewidth. Since the pressure broad- 
ening coefficient of individual absorption line for this gas sample is 



unknown, this comparison gives the correctness for the relative vari- 
ation between neighbouring lines rather than the correctness of the 
absolute value for absorbance. The relative variation is reasonably 
well reproduced and is limited by the signal-to-noise ratio of the 
measurements and the accuracy of the line strength calculations 
included in the database. 

Discussion 

The proposed method can enhance the spectral sampling interval to 
the level of the intrinsic linewidth of the comb mode (typically, less 
than a few Hz) 24 by increasing the observed time window of the THz 
pulse train and decreasing the incremental step of the comb sweep- 
ing, in principle. However, the actual linewidth in the comb-mode- 
resolved spectrum may be limited by the timing jitter between the 
two fibre lasers. Here we discuss the contribution of the timing jitter 
to the linewidth of the comb mode. This contribution is estimated by 
the accumulated timing jitter per time increment at the end of the 
sample interval and the number of the comb modes because such the 
timing jitter causes fluctuation in the size of the time window and 
hence the frequency increment in the Fourier spectrum 30 , which 
leads to broadening of the comb mode linewidth. In the case of the 
spectrally interleaved THz comb with a mode linewidth of 2.5 MHz, 
the temporal waveform of 100 consecutive THz pulses was acquired. 
The accumulated timing jitter was 117 fs for the 100-th THz pulse 
when the temporal waveform was measured at a sampling interval of 
25 fs; this contributes to a 1.75 kHz increase in the mode linewidth at 
0.6 THz. Therefore, we can conclude that the comb mode linewidth 
in the present system is not influenced by the timing jitter but is 
determined by the size of the time window; this is demonstrated by 
the absorption linewidth in the results shown in Fig. 4(d). This indi- 
cates that the present system has further linewidth gains that can be 
implemented. 
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Table 1 Spectroscopic performance of various THz spectrometers 
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Next, we discuss the applicability of the spectral interleaving tech- 
nique to an optical comb. In the case of a THz comb, the only way to 
achieve the spectrally interleaved comb is by sweeping f repl because 
the THz comb is a harmonic frequency comb o£f repl an d does not 
contain f ce0 . Sweeping off repl has the advantage that a small incre- 
ment oif repl causes a large change in the comb mode frequency due 
to the large number of comb modes. However, the resulting comb 
modes are not always interleaved appropriately over the full spectral 
range but are partially interleaved around the target frequency 
because the amount of frequency shift of each comb mode depends 
on its mode number m as well as Sf repl [see Fig. 1(c)]. For example, 
the comb gap is more densely interleaved below the target frequency 
whereas it is more coarsely interleaved over the target frequency. On 
the other hand, since the absolute frequency of the optical comb 
mode is given by/ ceo + mf repl , it can be spectrally interleaved by 
sweeping f ce0 as well asf repl . Although the previous studies of inter- 
leaved comb in optical region were based on sweeping of f repl 13 ~ 16 , 
sweeping off ceo allows us to make the optical comb interleave uni- 
formly over the full spectral range. One technical challenge of this 
method is large tuning off ceo overf repl /2 in order to fully cover the 
frequency gap between comb modes. Extension of the fully inter- 
leaved technique to an optical comb will further expand the possible 
applications of the optical-comb-based, broadband spectroscopy due 
to the uniformly spectral sampling density. 

Conclusion 

To achieve a spectral sampling density equal to the linewidth of each 
comb mode, frequency gaps between THz comb modes were success- 
fully interleaved by using swept dual THz combs. This is the first 
demonstration of overcoming the inherent limitation in THz comb, 
namely, the excessively discrete distribution of the comb modes 
limiting the fine spectral sampling for broadband spectroscopy. To 
demonstrate its superiority, we applied the interleaved THz comb to 
high-precision broadband spectroscopy of gas-phase molecules at low 
pressure. The resulting spectral resolution and accuracy were 
2.5 MHz and 8.39 X 10~ 7 , respectively. The proposed method has 
the potential to achieve the precise analysis of gas molecules in the 
atmospheric and space environments as well as combustion processes 
because the rotational transitions give them particularly rich spectral 
fingerprints in THz region while THz radiation is insensitive to scat- 
tering particles in optical region. Simultaneous achievement of high 
resolution, high accuracy, and broad spectrum by spectrally inter- 
leaved dual-THz-comb spectroscopy would enable us to discriminate 
densely distributed absorption lines correctly even though the target 
gases are mixed with aerosols, smoke, dust, clouds, or soot; for 
example, gas analysis in the smoke 31 or sooty flame 32 . 



Methods 

Experimental setup. The individual repetition frequencies of the dual mode-locked 
Er-doped fibre lasers (f repl = 250,000,000 Hz and/ rep2 = f repl + f oJfset = 
250,000,050 Hz) and the frequency offset between them (f offset = f rep2 ~ frepi = 
50 Hz) were stabilized by two independent laser control systems referenced to a 
rubidium frequency standard (Rb-FS, accuracy — 5 X 10" !1 , instability — 2 X 10~ n at 
1 s), as shown in Fig. 2. Furthermore, f repl and/ rcp2 could be respectively tuned over a 
frequency range of ± 0.8% by changing the cavity length. After wavelength conversion 
of the two laser beams by second-harmonic generation (SHG) with nonlinear crystals, 
pulsed THz radiation was emitted from a dipole- shaped, low- temperature -grown 
(LTG), GaAs photoconductive antenna (PCA1) triggered by pump light {X c — 
775 nm, At — 80 fs, P mean = 19 mW), passed thorough a low-pressure gas cell 
(length — 500 mm, diameter — 40 mm), and was then detected by another dipole- 
shaped LTG GaAs photoconductive antenna (PCA2) triggered by probe light (a c — 
775 nm, At — 80 fs, P meatl — 9 mW). The optical path, in which the THz beam 
propagated except for the part in the gas cell, was purged with dry nitrogen gas to 
avoid absorption by atmospheric moisture. A small portion of the output light from 
the two lasers were fed into a sum -frequency- generation cross-correlator (SFG-XC). 
The resulting SFG signal was used to generate a time origin signal in the ASOPS-THz- 
TDS. After amplification with a current preamplifier {AMP, bandwidth = 1 MHz, 
gain — 4 X 10 5 V/A), the temporal waveform of the output current from PCA2 was 
acquired with a digitizer (sampling rate — 2 X 10 5 samples/s, resolution — 20 bit) by 
using the SFG-XC's output as a trigger signal and the frequency standard's output as a 
clock signal. Then, the time scale of the observed signal was multiplied by a temporal 
magnification factor oi f repl /f offset ( = 250,000,000/50 - 5,000,000) 18 . This sampling 
rate and this temporal magnification factor enabled us to measure the temporal 
waveform of the pulsed THz electric field at a sampling interval of 100 fs. 

Comparison with other THz spectrometers. To underline the principal advantage of 
the spectrally interleaved dual-THz-comb spectroscopy, we compared its basic 
performance with that of conventional spectrometers using tunable CW-THz wave 
and broadband THz wave. Table 1 summarizes the comparison of the frequency 
coverage, frequency resolution, frequency accuracy, and dynamic range of power 
(DR) among them. CW-THz spectrometers using frequency- multiplied microwave 
sources (FMMS) can achieve the high spectral resolution and accuracy at high DR by 
referencing a microwave frequency standard whereas tuning range of a single FMMC 
is usually limited within 10 to 20% of a centre frequency 33,34 . CW-THz spectrometer 
with backward wave oscillator (BWO), also called to fast-scan submillimeter 
spectroscopy technique (FASSST), can extend the frequency coverage more widely 
while maintaining the similar resolution; however, its frequency accuracy is limited 
by the performance of Fabry-Perot cavity for frequency monitoring of BWO 35 . 
Although difference frequency mixing of C0 2 laser and 15 NH 3 laser with tunable 
microwave sideband generation (TuFIR) provides a spectral resolution below 20 kHz 
for CW-THz spectrometer, great care is required in obtaining precise frequency 
calibration and coherence 36 . Furthermore, in those CW-THz spectrometers above, it 
is difficult to fully cover the frequency range from 0. 1 THz to a few THz at once; it is 
necessary to select a suitable frequency multiplier chain, BWO tube, or laser emission 
frequency depending on the frequency band used, due to the limited tuning range 33 " 36 . 
Photomixing of tunable near-infrared lasers is a promising approach to achieve 
continuous tuning over 1 THz while maintaining frequency resolution of 1 MHz 37 . 
However, the frequency accuracy depends on the performance of a wavemeter that 
monitors the wavelength of near-infrared lasers. Although combination of the 
photomixing with optical comb makes the output frequency traceable to a microwave 
frequency standard, the continuous tuning range is largely reduced 38 . THz quantum 
cascade laser (THz-QCL) has the intrinsic characteristic of narrow linewidth 39 as well 
as high output power 40 . Furthermore, introduction of external cavity configuration 
provided broad tunability to this source 41 ; however, its continuously tuning range is 
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limited within 10 GHz whereas its discontinuous tuning range is achieved over 
100 GHz. On the other hand, phase-locking of THz-QCL to the optical comb enables 
us to obtain the frequency accuracy traceable to a microwave frequency standard and 
the reduced linewidth 42 . However, the continuous tuning range is limited within 
1 GHz 43 . 

In THz-TDS, the spectral resolution is determined by the length of mechanical 
time-delay scanning whereas the spectral accuracy depends on the positioning 
accuracy of a mechanical time-delay stage 44 . In ASOPS-THz-TDS, the spectral 
resolution depends on the laser repetition frequency: a lower repetition frequency 
leads to lower spectral resolution 20 whereas a higher repetition frequency contributes 
to higher DR and wider spectral range 19 . On the other hand, since the spectral 
resolution of the spectrally interleaved dual-THz-comb spectroscopy depends on the 
comb mode linewidth, better spectral resolution can be achieved regardless of the 
repetition frequency. The frequency accuracy in ASOPS-THz-TDS and interleaved 
dual-THz-comb spectroscopy is determined by the spectral sampling interval and the 
stability of a temporal magnification factor of f repl /j ^ 0 ff set , which is stabilized by the 
precise laser control. The high total balance of frequency resolution, frequency 
accuracy, and frequency coverage makes the interleaved dual-THz-comb spectro- 
scopy more attractive than other THz spectrometers though its DR should be 
improved to widely apply it to various spectroscopic analyses. 
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